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ABSTRACT 

The Stark-indnced shift and asymmetry, the so-called pressnre shift (PS) of 
Hq, and Balmer lines in spectra of DA white dwarfs (WDs), as masking effects 
in measnrements of the gravitational red shift in WDs, have been examined in 
detail. The results are compared with our earlier ones from before a quarter of 
a century (Grabowski et al. 1987, hereafter ApJ’87; Madej & Grabowski 1990). 

In these earlier papers, as a dominant constituent of the Balmer-line-prohles, 
the standard, symmetrical Stark line prohles, shifted as the whole by PS-effect, 
were applied to all spectrally active layers of the WD atmosphere. At present. 
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in each of the WD layers, the Stark-line-prohles (especially of H^) are imma- 
nently asymmetrical and shifted due to the effects of strong inhomogeneity of 
the perturbing helds in plasma. To calculate the Stark line-prohles in successive 
layers of the WD atmosphere we used the modihed Full Computer Simulation 
Method (mFCSM), able to take adequately into account the complexity of local 
elementary quantum processes in plasma. In the case of the Hq, line, the present 
value of Stark-induced shift of the synthetic Hq, line-prohle is about twice smaller 
than the previous one (ApJ’87) and it is negligible in comparison with the grav¬ 
itational red shift. In the case of the line, the present value of Stark-induced 
shift of the synthetic line-prohle is about twice larger than the previous one. 
The source of this extra shift is the asymmetry of peaks. 

Subject headings: physical data and processes: atomic processes, line: formation, 
line: prohles, plasma — stars: white dwarfs 


Introduction 


Quarter of a century has passed since we got to know for certain that the Stark-induced 
shift and asymmetry effects (hereafter called the “pressure shift”, PS), caused by local atomic 
quantum processes in plasmas, play an important role in formation of the Balmer hydrogen 
lines in white dwarfs (WDs) spectra, and that contribution of PS to the residual (free of the 
Doppler component) red shift of these lines must be taken into account while separating the 
gravitational (Ein s tein’Q red shift. We rn e et the hrst note on that theme in the papers by 
Wiese fc Kelleherl (1197111 and IWiese et ahl (1197211 . They experimentally revealed systematic 
slight shifts of these Stark-broadened lines to the red, which are linear functions of the 
electron density. As the source of light they used hydrogen plasma — the environment 
similar to that in the white-dwarf atmospheres. Moreover, the measu rements were taken in 
the sa me w ay as those in the WD s pectr a in two pioneering papers bv iGreenstein fc Trimble 
(119671) . and iTrimble fc GreensteinI (ll972l) . where the gravitational red shifts were determined 
within the relatively broad central region of a Balmer line — the so-called “averaged line 
center” (ALG) — in the best spectra available at that time, i.e. in those recorded by means 
of a 508-cm telescope. 

The results of the gravitational red shifts determination by Greenstein and Trimble for 
white dwarfs r evealed an unexpe cted and surprising problem, with time called the astrophys- 
ical puzzle (cf. IWeidemanrull979|) : due to some unknown reasons, the measured (residual, i.e. 
free from the Doppler component) WD red shifts were systematically too great, up to about 
50 per cent (10 - 15 km/s in Doppler velocity units) exceeding predictions by Einstein’s 
















3 


General Relativity (GR). (More precisely: WD masses resulting from the gravitational red 
shift measurements systematically exceeded those determined usin g the so-called astrqphys i- 
cal methods, i.e. astrometric one, model atmospheres, and so on.) IWiese &: Kelleherl fll97ll) . 
in their paper under the significant tittle “On the cause of the redshifts in white-dwarf spec¬ 
tra”, noted a propos: “It appears that the Stark-induced red shift accounts for a portion 
of the observed red shift in white dwarfs which has previously been attributed entirely to 
gravitation.” 

At that time, the appearance of the Stark-induced red shift was surprising indeed, 
because splitting of Balmer lines in external homogeneous electric helds, as observed in 
those days in laboratories, was symmetrical. Quadratic Stark e ffect generates line shifts 
(usually to the red) and line asymmetry. “However — we quote IWiese &: Kelleherl fjl97l[) 
ag ain — based on all the then available experim ent al and t h eoret ical knowledge — the works 
of iRausch von Traubenberg fc Gebauerl (119291 ) and iGriemI (119641) are especially cited — the 
Stark shifts were estimated to be insignificant.” 


In both objects under consideration — in the WD spectra and in the light of the 
device (the high-current, wall-stabilized arc) of Wiese-Kelleher’s (1971) type — the hydrogen 
Balmer lines are enormously broadened by the interparticle Stark effect in plasma. In both 
cases also these lines, for a long time, appeared to be practically symmetrical and unshifted, 
contrary to the so-called “metal” lines of all heavier atoms — narrow and easy to measure, 
but strongly affected by Stark shift and asymmetry. These last (“metal”) lines are completely 
useless for the purpose of the gravitational red shift measurements. For about four decades 
now it has undoubtedly been a well-known fact that the Stark-broadened hydrogen Balmer 
lines in dense plasmau are — apart from their enormously great broadening — slightly 
asymmetrical and red-shifted as well. The asymmetry and red shift are, relatively, of small- 
size effects, but phenomenally, in their physical specification, very complicated ones. Detailed 
measurements and a careful theoretical balance of physical processes forming the hydrogen 
spectral lines in plasma step by step lead us to complete understanding and to a successful 
description of the Stark-induced shift. So, the Stark-induced shift, as a masking effect, can 
be definitely separated from the gravitational red shift measurements in the WD spectra. 


The first simple attempt to provide a theoretical explanation of the effect discov¬ 
ered by Wiese et al. (1971; 1972) in the hydrogen Balmer lines in plasma was made by 


^One routinely reads the term “dense plasma” (or “high-density plasma”) as a plas ma of electron con- 


Wiese et all (Il972h . entitled “Detailed 


centration about and above 10^® cm“®. For example, in the paper by 
Study of the Stark Broadening of Balmer Lines in a High-Density Plasma”, investigations cover the range 
of electron densities W between 1.5 x 10^® and 10^^ cm“®. In WD atmospheres, the range of Ne is usually 
considerably broader. 
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Grabowski &: Halenka fjl975h — as a specific consequence of Debye screening, i.e. of the 
neutralizing action of a negatively charged plasma background on the atomic nucleus, and 
of resulting deformation of the internal structure of disturbed hydrogen atoms. The result, 
PS oc iVe/T (for details see that reference), was surprisingly fair: the predicted pressure shift 
of the line center to the red became evident and was of the same order of magnitude as the 


observ ed excess over gravitational red shift. In the recapitulating paper bv iGrabowski et ah 


(1l987l) . the synthetic Balmer line-profiles in WD spectra were calculated for a wide range 
of log^f and values of WD atmospheres in a physically more realistic approach, cor¬ 
responding to the world state of the art in the hydrogen-spectral-line broadening theory 
of that time (see Section p). W e used, namely, symmetrical Stark-broadening-line-profiles 


ffollowing lEdmonds et al.lll967() as functions of the local physical conditions in subsequent 
layers of a stellar atmosphere, in each layer, moreover, shifted as a whole on the value of 
the pressure shift, PS, dependently on proper local physical condition in the layer. The 
PS-function was calculated, taking into account complex influences of plasma environment 
on a H atom: (1) predomi nating effect — the screening o f the atomic nucleus by the cloud 


of free electrons, following iGrabowski fc Halenka (1197511 . but corrected empirically (by a 
multiplicative factor), accordingly to the laboratory data accessible at that time, and (2) the 
phenomena of the correcting nature (especially in the case of the line), caused by con¬ 
tribution of the resultant quasi-static local ionic held. We took into consideration following 
contributions: multipole expansion of the perturbing helds, the quadrupole term including, 
and the second order of ionic held’s corrections in frequency shift (quadratic Stark effect). 
The hrst order correction in the oscillator strength alterations of the line components, and 
the so-called pres s ure io nization in plasma, have been taken into account also (for details see 
Grabowski et ah ( 1987 )1. 


The principal results were as follows: Asymmetry and pressure shift (to the red) of the 
synthetic Balmer line increase with the distance from the line center; the line is strongly 
affected by PS and, consequently, is less useful for the gravitational red shift determination, 
especially when far wings are used to dehne the averaged line center. The relative signihcance 
of the pressure shift of Balmer lines, as compared with the Einstein’s gravitational red shift, 
systematically decreases with increasing values of Teg and of logg — from about 50 per cent 
of the gravitational red shift down to about mere 1 percent -— at both opposite poles, i.e. 
at extremely low and extremely high values, respectively (cf. iMadei fc Grabowski! (119901) ). 


Today, over a quarter of a century later, when excellent, high-resolution spectroscopic 


HST 


T.g. 


Barstow et ah 

2005 

), ESO (e.g. 

Falcon et ah 

2010 

), or VLT (e.g. 

Gasewell et ah 


20091) . the importance and usefulness of these old, above-discussed, results on PS are rather 
archival. At present we are able to calculate the synthetic WD Balmer line-profiles with 
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incomparably higher precision, taking into account a more complete description of the quan¬ 
tum physical processes, active in the WD Balmer line creation. (For comment on the specific 
phenomena, active in creation of PS and asymmetry of Balmer lines in plasma, taken into 
account in the present calulations, see Sec. 2, the second through the last paragraphs, 
especially the visualization in Figs. [1] and [2l and also in Sec. [3l the third and the fourth 
paragraphs.) 


Measurements of the gravitational red shift give WD masses (precisely: mass/radius 
ratio) immediately and with the highest exactness (of the same order as — in the incidental 
cases — the orbit astrometry of binary systems), practically independently of anything 
— under condition, however, that the spectral measurements are free (separated from) of 
disturbing and masking effects, produced by mixed up elementary processes in plasma of WD 
atmospheres, especially at low values of log^f and Tgff. We devote our paper to improving 
the relia bility of this excellent method. Activity in that direction is desired and exp ected 


(cf. e.g. iKoesten 119871: iProvencal et al.l Il998l: iBarstow et al.l l2005l: iFalcon et al.l l2010l) . We 


will discuss relevant problems in next sections. 


2. An introductory note on the physical elementary processes forming the 
Balmer line profiles in plasma (in WD atmospheres in particular) 

One cubic centimeter in a sample layer of the typical WD atmosphere, at temperature 
of the order of ten thousand Kelvin, contains about 10^° atoms, almost all are hydrogen ones, 
of which one in a thousand is an ionized atom, i.e. without a bound electron. So, in that 
volume we meet about 10^^ electric charges (electrons and ions). The electric field, attached 
to a charged particle, is a long range one (~ r“^), therefore each charged particle acts 
simultaneously on lots of other charges and neutrals. All particles run with great velocities, 
of hundreds/thousands of meters per a second. In the scale of nanoseconds (the typical clock 
beat for atomic micro processes), the influences by ions (slow), are nearly constant, resulting 
in quasi-static splitting of the atomic levels into sets of sublevels; the perturbations by 
electrons (fast) — rapidly changing, are “seen” by atoms as impacts, resulting in a reduction 
of life-times and in broadening of the atomic sublevels. Moreover, mutually overlapping field 
strengths reduce the scope of this field penetration (the so-called screening effects). This is 
a typical scenario of line profile formation in the plasma conditions. 

For start, we consider the simplest, virtual case: a single ion (proton) at rest, acting on 
a tentative hydrogen atom, also fixed, dependently on the distance between these particles. 
We temporarily ignore all the remaining charged particles. This pair of particles (called the 
nearest neighbor approximation, NNA) stands for an extremely simplihed physical quantum 
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system, but the mathematical problem of so lution of the corresp onding time-independent 
perturbation theory is by no means simple fcf. lOlchawa et al.ll200ll. and Figs.[T]and|2]therein, 
and below). 


Fig. [T] presents the solution for the Balmer-a line as an example. In the homogeneous 
external electric held (in our demonstrative model of the perturbing point charge, the ho¬ 
mogeneous held formally corresponds to the inhnite distance of the perturber) the splitting 
of line components is strictly symmetrical (see up-plane of the hgure). This is the canonical 
pattern of the linear Stark ehect (LSE) of Balmer-a in a homogeneous external held. The 
held of the perturbing ion is, of course, the more inhomogeneous, the nearer the ion resides 
(the ion-atom distance is expressed in R/ao, oq — Bohr radius) and, consequently, the mul¬ 
tipole expansion (ME, some kind of the Taylor series) of the perturbing held strength must 
be taken into account in the calculations, i.e. the dipole term (corresponding to the ho¬ 
mogeneous component), a quadrupole one (the hrst order correction), an octupole one (the 
second order correction), and so on, if necessary. In the case under examination, all these 
higher order multipole terms have been taken into account. We see the ehect of these terms 
on the line splitting in the changes in location (and even in succession!) and in intensities 
of the line components. Particularly, the intensities of the red components (the right hand 
side of the hgure, when we are looking toward the LSE pattern) decrease with the decreasing 
ion-atom distance, but all line components move to the red — creating a complicated image 
of the pressure asymmetry and shift formation. 


In Fig. m the evolution of the intensities of the selected Stark components of the 
line is presented demonstratively. As a sample, the strongest (vrg) Stark components are 
used: 7r_8 (the blue component, the upper curve of the given mark in the hgure), and tt+s 
( the red component, the lower curve in the hgure), in dependence on the approximation 
used, as a function of the reciprocal reduced distance {ao/R). The marks used have the 
following meaning: solid lines — dipolar perturbation (homogeneous external held), the 
quadratic Stark ehect is included; the short-dashed curves — the quadrupole interaction is 
also included; the dotted curves — the octupole interaction is also included. The evolution 
presented by the dotted curves (marked by proper arrows) is the most precise. 

In Fig. [T] and, especially. Fig. |21 we see that the main source of asymmetry — and of 
the asymmetry-induced shift — in the quasi-static approximation is the inhomogeneity of 
the ion microheld in the surroundings of an atom, active in the spectral line formation: the 
quadrupole and the octupole terms, and also the quadratic Stark ehect. 





7 


3. Calculations of the Stark profiles of hydrogen lines in plasma environment 


The fundamental principles of the Stark-prohle calculations of hydrogen lines, formed in 
the plasrna envi ronments, were formulated in the 1950s, in the pioneering pape rs by Michel 


Barangerl (119581) and in the 1960s by Hans R. Griem and by his coworkers flGriem et ah 


19621) . In the so-called generalized impact-broadening theory, formulated at that time, the 


influences of the fast moving electrons perturb, as impacts, a hydrogen atom and reduce 
life-time (causing broadening) of sublevels (separated and well-dehned in quasi-static ion 
held) of both the upper and the lower levels of the line. A complete Stark prohle is a result 
of an averaging of impact prohles created in such a manner over the quasi-static ion held 
strengths within the Debye’s sphere (some “interacting volume” of plasma round an atom, 
limited by ehects of mutual screening o f the charged particle s). The best sources of the 
Stark prohles of those days were these bv iGrieiru fll964l) . 
fjl968l) . generall y known as “ modih ed impact theory' 


GriemI fjl974l) , and iKepple &: Griem 
Nearly equiv alent ones are: the “un i- 


hed theory” by IVidal et al.l (119731) , and “semiempirical prohles” by Edmonds et al.l (119671) . 
They account for the Debye’s shielding ehects of the perturbing ions, but the ions everywhere 
are treated as static ones (the thermic ion dynamics is omitted generally), and, furthermore, 
the Stark splitting caused by ions is assumed as LSE entirely — a symmetrical and unshifted 
one — i.e. the higher order corrections in the multipole expansion, ME, and in the pertur¬ 
bation theory, PT (the second order PT-correction gives the quadratic Stark ehect, QSE), 
are completely neglected. In all these calculations, the classical path approximation, GPA 
(assuming that the perturbing electrons move along the classically specihe d ways) is al ways 


used. Fortunately, GPA is entirely satisfactory in most applications (e.g. lGrieirull974j) and 


does not demand improvement. Negligence of the ion dynamics creates, however, serious un¬ 
certainty in the line centers — in the region of special importance for gravitational red shift 
measurements. An extreme defect — from the point of view of the gravitational red shift 
measurements in WD spectra — of the quasi-static approaches in the majority of existing 
Stark-prohle theories lies in the complete negligence of the higher order corrections in the 
Hamiltonian (ME) and in PT, as well. It is these corrections that produce the red shift and 
asymmetry features in the hydrogen-line prohles, observed in spectra of moderately dense 
and in high-density plasmas, and, naturally, appearing in the WD spectra, too. 


The spectral line shapes, particularly the hydrogen line shapes, are one of the most 
important carriers of information on the structure of laboratory, as well as on space plas¬ 
mas. From the astrophysical point of view, it refers mainly to the plasma of the stellar 
atmospheres and the plasma of WD atmospheres. In order to get reliable recognition of the 
physical features of the investigated plasma environment (i.e. in order to get the plasma 
diagnostics) through the line-shape analysis, the line shapes must be calculated with suffi¬ 
cient accuracy. To attain this different approaches are accepted — analytical methods and 































computer codes of various complexity and accuracy, and of different limits of applicability. 
The reader, interested in this subject, can find more informati o n on the line broadening in 
recent reviews, e.g. bv lGiinter &: KoniesI fjl997|l . iDiurovic et ahl (1200911 and in papers quoted 
therein. Particularly, one can hnd a review of the line prohles obtained by using the so-called 
Full Computer Simulation Method (FCMS) compared with other approaches in the paper 
by Ferri et ah (2014). 

The line shape and/or the line width (full-width-at-half-maximum, FWHM, in par¬ 
ticular) serves the purpose of plasma diagnostics, as t he main informer on the p h ysical 


properties of a diagnosed medium. FCSM approaches fICigosos &: Cardenosol Il989 


1996 


Halenka &: Olchawal Il996h very faithfully reconstruct the local phenomena in plasma and, 
consequently, give shape and width of the line relevant to this aim. These approaches are, 
however, able to predict neither asymmetry nor a shift of the line. 

We refer shortly to the principle of FCSMs. The local ionic, as well as electronic, shares 
in the spectral line formation are appointed using the simulation method, and the time de¬ 
pendent Schrodinger equation is solved numerically. The ionic contribution is, as a rule, 
limited to dipole interaction of the atom-ionic field, and to the first order correction in FT, 
i.e. to the linear Stark effect (LSE) only. The resulting line-prohle is there fore a symmetrical 
and un shifte d one. However, th e FCSMs have no limitations present, e.g., in iKepple &: Criem 
(119681) . and IVidal et ahl (119731) theories, i.e. quasi-static ion helds, nor an impact approx¬ 
imation for free electrons. The quadrupole interactions (the hrst order correction in the 
multipole expansion, ME), as well as the quadratic Stark effect (QSE, the second or d er cor - 
rection in FT) can be taken into account in FCSM, indeed — see paper bv I Olchawal (120021) . 
Such an improved approach leads to effects of the line-asymmetry and to Stark-induced shift 
- to effects important especially in conditions of dense plasma. These predictions excellently 
agree with the recent laboratory measurements. Unfortunately, this improved approach is 
very time-consuming, and a laborious one. We, therefore, have developed an equally reliable, 
but more effective method, useful in efficient line-profile calculations, when a wide range of 
different physical conditions (successive layers of the WD atmosphere) must be examined. 
Here the shape of the spectral line is a sum of two components. The hrst component (of 
the zeroth order term in FT and the hrst in ME, i.e. LSE in dipole approximation), is 
calculated using the FCSM approach. The second component (consisting of a variety of 
corrections, each small compared to the hrst component) includes the contributions origi¬ 
nating from higher-order terms of the emitter-ionic microheld interactions (the quadrupole 
interactions and QSE, calculated in quasi-static approximation for ions). This second com¬ 
ponent describes the shift and asymmetry of the Balmer line prohle, produced by the ion 
constituent o f plasm a . The shifts caused by electrons have been taken into account, following 
the paper bv ICriemI (119831) . The method used - modihed FCSM (in our notes, mFCSM) - 
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is an effective one, producing realistic line-profiles: broadened, shifted, and asymmetrical, 
agreeing, moreover, wi th the laboratory m e asurements. T his method has been formulated 


in a series of papers b y IHalenka fc Olchawal fll996l) , lOlchawal (120021) , IWuiec et ahl (120021) and 
Olchawa et al.l (120041) . and, hnally, was successfully examined and improved in the paper by 
Griem et ah ( 20051) . All the calculations of the hydrogen line-profiles in the present paper 
have been carried out following this mFCSM method. 

Particular problems, related to the matter in the hand, have been solved in the following 
references: A detailed discussion and proper quadrupole and octupole corrections of the ionic 
helds, producing asymmetry and red shift of the Stark line profile in plasma — in a series of 
papers by Halenka and Olchawa: IHalenka! (jl990l) . where correction of the “modihed impact 
theory” was made, including the er nitter-ion-quadrupole intera ctions, the quadratic-Stark- 
effect, and the ion dynamics effects; iHalenka fc Olchawal (l2007l) — octupole inhomogeneity 
tensor o f ionic microfield in Debye plasma at neutral emitter, calculated for the hrst time; 
Halenka! (120091 1 — the octupole inhomogeneity tensor at ionized emitter, also calculated for 
the first time. 


The starting point of this consideration is the relation between the spectral line profile 
and the average of the dipole autocorrelation function C{t), which can be written in the 
following way: 

I{Au)= lim Ti~^ [ C'(t)e‘^‘^Mt (1) 

Jo 

C{t) = Tr{di/ • d/v Uifi{t)}av/ Tr{dj/ ■ d/J, (2) 

where d is the dipole operator for the hydrogen atom, while ii' and ff indicate the sublevels 
of the initial {Ei) and final {Ef ) states of the unperturbed atom, respectively. The frequency 
separation from the line center is given by Au = u — {Ei — Ef)/h, whereas U{t) is the operator 
of the time development of the hydrogen atom in the presence of the electric held produced 
by electrons and ions. The averaging {• • -^av is taken over all initial simulated held strengths 
and possible time histories. The time-evolution operators Ui/iit) and Ujj:,{t) (corresponding 
to the initial and hnal states, respectively) satisfy the following Schrodinger equation: 

thU{t) = [Ho + V{t)]U{t), (3) 


where Hq is the Hamiltonian of the isolated radiator and V (t) is the radiator-plasma inter¬ 
action potential. The Hamiltonian, with th e multipole exp ansion of the potential restricted 


to quadrupole terms, may be written, e.g. 


3 multipole exp a 
Halenkal (jl99ol) . 


as follows: 


H{t) = Ho-d- F{t) - 1/6 + • • •, 

jk 


( 4 ) 
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where the second term in Eq. 0] describes the plasma-emitter dipole interaction and the 
third represents the quadrnpole interaction of the radiator with the inhomogeneons electric 
microheld of the plasma. The electric held F(f) and the inhomogeneons microheld tensor 
Fjk{t) represent the respective total helds originating from ions and electrons. 


We computed a pure H model atmosphere of typical WD parameters: Teg = 12 000 
K (in hgures below we use the abbreviation 12 kK) and \ogg = 8.0 (cgs units) using the 
ATM computer code. Synthetic Balmer line prohles were computed using our LINE code, 



present paper. Numerical accuracy of the ATM code is high, since it ensures that error 
of the bolometric radiation hux is much less than 1% across the WD model atmosphere. 
Both high accuracy and the physical foundations of these codes (requirement of hydro¬ 
static and radiative equilibrium) are comparable with other model atmosphere c omputer 
programs used routinely now adays, e.g. the computer code TLUSTY Version 195 (iHubenv 


19881 : IHubenv fc Tja.n7Jll997|l in its basic component, modelling plane-parallel, horizontally 


homogeneous stellar atmosphere in radiative and hydrostatic equilibrium. LINE, which is a 
program for calculating the spectrum emergent from ATM m odel atmosphere, is comparable 
with program SYNSPEC Version 43 flHubenv fc LanzIPfinnh . The only serious simplihcation 
is the LTE solution of the equation of state. Therefore, both ATM and LINE codes are 
suitable for differential investigations of the features of the synthetic Balmer-line profiles, 
when our interest is focused on testing of two Stark opacity functions, distinctly differing 
from each other. 


Our model atmosphere of the sample white dwarf consists of 49 gas layers of various 
optical and geometrical depth. The coolest (shallowest) layer has temperature T = 5 500 K 
and electron concentration about W = 1.5 x 10^° cm“^, whereas the deepest (hottest) 
one, contributing to the synthetic hydrogen spectral line prohle, has T = 55 000 K and 
Nf, = 2 X 10^® cm“^. For each layer of the WD model atmosphere we calculated complete 
Stark line opacity profiles, which were broadened, asymmetrical, and shifted as the result 
of all relevant physical processes in plasma, including motion and shielding of perturbing 
ions plus the Stark-Doppler coupling (not only convolution, but also dependence). For this 
purpose we applied the modified FCSM. 




















11 


4. Verifying the theoretical calculations 


In order to check the calculations of the complete — broadened, asymmetrical, and 
shifted — Stark prohles of the Balmer-a and Balmer-/3 lines in the proper range of the 
physical conditions of hxed, homogeneous, thin-layer plasmas, we have carried out verify¬ 
ing experimental measurements. For the purpose of this study we have applied the same 
experimental device of the plasma generation (a high-current wall-stabilized arc), the light 


treatment, and the smoothing p roced ure of the spe c 


her papers: IWuiec et al.l (1200211 . and iGriem et ahl (120051) . Srj ecihc data of the experiments 


were nearly the same as in the paper by iGriem et ah 


ral re cordings, as those described in ear- 


(1200511 


with differences concern¬ 
ing the working gas (argon-hydrogen mixture in the cited reference) and the method of the 
plasma diagnostics (LTE), as follows: (1) As a plasma material we have used nearly pure 
argon, with a small admixture of the atomic hydrogen gas. (2) The pLTE (partial LTE) 
plasma diagnostics have been applied. 

In the cylindrical channel of the discharge arc, working at the atmospheric pressure, 
the argon plasma is very stable, of nearly perfect cylindrical symmetry, and homogeneous 
along the lines of view parallel to the arc axis. The argon plasma, with a small admixture 
of hydrogen, acquires features of high-density plasma, and, at the same time, the plasma 
is an optically thin layer in hydrogen lines. T and N^. values are extremely high along the 
symm etry axis of the arc channel, and decrease with the radial distance from the axis. In the 
paper IGriem et al.l (120051) . we diagnosed plasma conditions assuming local thermodynamic 
equilibrium (LTE). In the present paper, however, we carried out — more careful and real¬ 
istic — pLTE plasma diagnostics, i.e. we adopted: modified Boltzmann’s and Saha-Eggert’s 
laws, Dalton’s law and the electrical neutrality of the plasma. We took the integral emission 
coefficients of Arl 4300 A and or from the experiment. We have take n the overpop¬ 
ulatio n factors of Arl ground level for pLTE diagnostics from the paper by iHelbig fc Nick 
(1l98ll) . The experimental investigations have been performed at plasma electron densities 
between 2.8 and 10.0 x 10^® cm“^. Some results of the experiments are visible in Figs. 3 - 
6, as an example. In the matter of pressure shift and asymmetry of the investigated Balmer 
lines, we obtained an excellent agreement of the experimental and the calculated data. 

In Fig. [3l the Stark line-profiles of the are compared. The black points — the exper¬ 
imental data of the present paper; the long-dashed curve — the same theoret ical approxima¬ 
tion a s that used in our paper ApJ’87 — symmetrical Stark profile following lEdmonds_et_aL 


(119671) . hereafter ESW, but shifted as a whole, in dependence on the physical conditions of 
pl asma in the g i ven la yer of the WD atmosphere; the short-dashed curve — Stark line-profile 
bv IVidal et al.l (1l973l) . hereafter VCS — a symmetrical an d unshifted one ; this approxima¬ 


tion is the main one in the astrophysical spectroscopy (e.g. iBarstow et al.l (1200511 ): the solid 
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curve — our modified-FCSM, complete Stark line-profile of Balmer Hq, asymmetrical and 
shifted to the red. We see particularly that: (1) Agreement of our calculated (using mFCSM) 
and measured Stark line-prohle of Balmer Hq, is a perfect one; (2) ESW — shifted (ApJ’87) 
— the line prohle runs correctly; (3) VCS line-prohle has a defectively described line center, 
and is not suitable for WD applications. 


The FWHM of the line and the red shift of the “experimental line center” (ELC, as 
dehned bv lWiese fc KellelieiillQTlh of the experimental line-prohle and of the modihed-FCSM 
one are also shown in the hgure. The agreement of these data in the case of the line, 
as well as the line (Fig. H] below), is excellent. We would like to pay our attention to 
the pressure shift of the Hq, line peak — a little smaller one in the experimental and in 
modihed-FCSM line-prohles than that in paper ApJ’87, similarly as in the case of the 
line, discussed in detail in Figs. E] and [H] below. 


In Fig. [5l we see that in ApJ’87 the PS-values of the peak of the Hq, Stark line-prohle 
are considerably overestimated. In the case of the line, the circumstances are diherent. 
In Fig. [HI it is clearly seen that PS of the central dip of Stark line-prohles were predicted 
well in ApJ’87, with a little excess as compared with the present, more precise data. The 
PS measured bv IWiese et al.l (119721) . presented in Fig. [6l require a comment. We see that 
their best-htting curve achieves zero electron concentration (i.e. the conditions of the iso¬ 
lated atom) at the blue pressure shift of about 0.2 A, which is, of course, not a reasonable 
result. Probably all these measurements of the line were taken using temporarily de¬ 
fected recording equipment, producing unphysical, systematic (additive) red shift. (We see 
that Wiese’s and our curves of the PS are nearly parallel.) When the presumed additive 
contribution (probably erroneous) is removed, the measurements by Wiese at ah precisely 
agree with our measurements and the modihed-FCSM calculations. 


5. The emerging (synthetic) WD line-profile calculations 

Complete — broadened, asymmetrical, and shifted — Stark line-prohles have been 
calculated for each layer of the WD atmosphere under examination, and, subsequently, 
convolved with the proper Doppler contribution. In Figs [7] and [H we see, as an example, the 
hnal convolution of Stark-Doppler ehect for one of the shallowest layers (a relatively cool 
and rare one) — Fig. [TJ and for one of the deepest (hot and dense) atmospheric layers — 
Fig. [HI 

The emerging (synthetic) Balmer line-prohle is a weighted mean, balanced by the trans¬ 
fer equation, of the contributions of all active, outer layers of the pure H model atmosphere. 
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The superposition of the complete Stark-Doppler prohles of the Balmer line (asymmet¬ 
rical and red shifted), corresponding to successive atmospheric layers — from the extremely 
outer layer (narrow structure) down to the deepest one (extremely broad shapes) — is shown 
in Figs. OandUnj (The H-alpha line was treated just in the same maner.) 


In Fig. ITTh. the synthetic spectrum in the region of the line, calculated using opacity 
prohles, as in Figs. |9] and [101 is presented in the huxes scale. This “wide-angle” spectroscopic 
view is, of course, very sir nilar to those observe d via most advanced pr esent-day instrumenta - 
tion, as reporte d, e.g., bylBarstow et ah! (120011), Fig. 01 ground-based; iBarstow et ah! (1200511 . 
Fig. 5, HST; or Kepler et ah (2007), Figs. 10 and 11, SDSS. In Fig. ITTb — the same, but in 
a reduced (0-1) scale, r\ = Fluxiine/Fluxcontinuum- From the point of view of the gravitational 
red shift measurements, the clue lies in the line center — in the narrow central section, barely 
a few Angstroms wide. We discuss the hnal results of the analysis of that region in detail in 
the next section. 


6. Results and conclusions 


The pressure shift in the central parts of the synthetic Ha Balmer line is practically 
negligible; our earlier (ApJ’87) diagnosis in that mater is hereby conhrmed. (N.B.: in 
very-high-density plasma conditions that scenery is more complicated — cf. the paper by 
Griem et al.ll2005l . e.g. Figs. 9 and 13 therein — the pressure shift of the line achieves 
substantial values.) In the case of the synthetic line, the central region is slightly affected 
by the pressure shift, as correctly predicted in ApJ’87, but the near and far wings are affected 
enormously, as compared with ApJ’87. It is seen in our Figs. 02] and [121 plotted on the basis 
of the synthetic spectra of the WD atmospheres of the typical parameters: log^f = 8.0, and 
Teff = 10 000 K (in hgures we have used abbreviation: 10 kK), as an example. 

We consequently use the concept of the pressure red shift of the WD Balmer lines as the 
red shift of the “averaged line center” (ALC), introduced bv [Greenstein &: Trimbld f|l967[) 
to the WD spectra, and ad apted in the exper i ment al measurements, as the ELC paradigm, 
by Wiese fc Kellehej f 197ll) and Wiese et ah ( 1972 ). So, in Figs. 02] and 031 on the abscissa 
axis the quantity AA (A) is used; on the ordinate axis — the pressure red shifts in the 
equivalent Doppler velocity units [vp or Vr, respectively), in conformity with the tradition in 
that subject. In Fig. 021 the lower index F at the velocity symbol reminds that the ordinate 
axis of the synthetic WD spectrum is the fluxes scale F\ (not the residual intensity r\). The 
thick curve in Fig. 021 has been constructed by us from the WD synthetic spectrum, strictly 
in the same manner as an observer does on raw recording of the real WD spectrum. Here 
we must note, however, that the straightforward measurement of the WD gravitational red 
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shift in ALC-convention can be falsified by apparent, “trivial” causes, the prin cipal of whic h 


is the red shift induced by the sloping of the Paschen continuum (cf., e.g., ISchulzl 119771) . 


The measurements in the residual intensities, r^, in the (0-1) scale, are free of that “trivial” 
contribution. 


In Fig.[T3l the main results of this paper — the Stark-induced red shift of the investigated 
Balmer lines for a WD of the given parameters — in the reduced scale, r;,, are presented. We 
see here four curves. One pair of the curves correspond to Hq,, another — to the line. The 
solid lines show new effects — the results of the present paper; the dashed ones refer to the 
results of ApJ’87, with the aim of comparing these results. Each of the shown dependences 
presents the Stark-induced red shift (pressure shift, PS) alone, i.e. net pressure line shift 
(free of “trivial” contributions), produced by physical, elementary processes in plasma of the 
WD atmosphere, in all its layers active in the formation — via the transfer equation — of 
the synthetic shape of a Balmer line. 

We remember here the fundamental differences in the approaches applied in the earlier 
papers ApJ’87, A&A’90 and in the present paper: In the former ones, the symmetrical 
line-prohles of the Balmer lines, shifted as a whole, were used. The pressure shift of the 
line center, used at that time, remains nearly correct. However, the synthetic line-profiles 
— calculated at that time using symmetrical line-profiles, shifted as a whole — have been 
uncertain. In the present paper, the Stark-line prohles used are not only red shifted but also 
asymmetrical ones, similarly as in real conditions. The results presented in Fig. ITS] permit 
us to formulate the following conclusions: (1) The Stark-induced red shift (pressure shift, 
PS) of the line in WD spectra must be taken into account when the gravitational red 
shift measuring, especially in cases when the ALC paradigm is applied. In the line in 
WD spectra, the Stark-induced red shift is far greater than that in the case of the line; 
in the case of the last one, the pressure red shift can practically be neglected. (2) Our new, 
values of the “observed” pressure red shift of the WD line are nearly twice smaller than 
those predicted in ApJ’87, whereas in the case of the WD line — considerably greater, 
especially at AA > 10 A, compared with those in ApJ’87. 

The effects of such clear differences in the observed features of both these lines, and 
in the resulting Stark-induced red shift in WD spectra, can be explained as a result of 
dissimilarity of the quantum structures of both lines. The line belongs to the so-called 
“odd” Balmer lines, with a strong central line-component, whereas the H/j line — to “even” 
ones, without the central component. These differences strongly manifest themselves in the 
laboratory (emissive, and thin-layer plasma; cf., e.g.. Figs. [2]and|l]) plasma spectroscopy. In 
the inhnitely thick stellar atmospheres, the “odd-even” problem is — from the point of view 
of the observer — completely masked by processes of the transfer of radiation to outer layers 
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of the stellar atmosphere and to the observer’s eye: the synthetic lines of both these kinds 
have similar shapes, with a clearly formed central dip. In each atmospheric layer, however, 
the line-shapes of opacities have snch featnres as in a thin layer — cf. Figs. 7-10 of the 
present paper. 

The H,g line-shape is a specihc one (see Figs. H] and |8]) — the red peak of the line is 
lower than the bine one (i.e., the “bine” opacity is greater than that in the “red” line region). 
Conseqnently, a radiation transfer in the region of red peak proceeds more easily compared 
with that in the bine peak. In the red side of the line, we “see” therefore more deep, dense, 
and hot layers; on the bine side of the line — more shallow, relatively rare and cool layers. 
These effects averaged (by a radiative transfer eqnation) prodnce additional (as compared 
with contribntion of a single atmospheric layer) asymmetry and Stark-indnced red shift of 
the synthetic line-prohle. The case of the Hq, line is, dne to a strong central component and 
nearly symmetrical Stark splitting, incomparably simpler; asymmetry effects (appearing in 
this line too), do not prodnce any snbstantial Stark-indnced red shift. 


When spectroscopic data of extremely high dispersion are available, e.g. from HST or 
VLT, and when the line with the central non-LTE dip (cf., e.g.. Falcon et ah 2010, Fig. 2, 
or Casewell et ah 2009, Figs. 2 and 3) measnrement is nsed, the problem of the Stark-indnced 
shift can be completely ignored. The non-LTE dip is formed in the nppermost, extremely rare 
atmospheric layers of WD, where Stark effects are of no importance as against Doppler effect 
(cf. Fig. [7]). Here the observed line-shift is entirely cansed by the gravitational (Einstein’s) 
effect. In more freqnent cases of the midd l e-resolntion spectr a, for example from SDSS, (e.g.. 


Kleinman et ahl 120041: iKepler et al.l 120071: iMadei et ah 


200411 . at the gravitational red shift 


measnrements the “average line center” dehned in broader regions of the line shonld be 
nsed. The same concerns the older archival spectroscopic data, often invalnable records of 
the past days. 


When the synthetic or Balmer lines-prohles in the residnal intensity nnits, r^, 
in the (0-1) scale are nsed (as we see, e.g., in Fig. [TTbl. the data from Fig. [131 free of all 
“trivial” agents and most reliable ones, shonld be nsed thronghont. 

We hope that the presented resnlts will be nsefnl in a more precise determination of the 
gravitational red shifts in the WD spectra and in a more accnrate determination of the WD 
mass/radius ratio. It is important from the point of view of onr nnderstanding of WDs — 
objects largely widespread in Natnre. Better knowledge of that ratio can also be helpfnl in 
imposing new limitations on the WD interior theory and can faciliate better calibration of 
la type snpernovas — exploding WDs in binary systems — and, conseqnently, can permit 
to more precisely examine the density of the dark matter in the Universe. 
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Fig. 1.— The quasi-static Stark spectrogram (pattern of splitting) of the Hq Balmer line, 
as an example, in the held of a single, hxed ion, as a function of the ion-atom distance R, 
expressed in the reduced scale R/ao, where Oq is the Bohr radius; at (formally) R/ao = oo 
(in the homogeneous perturbing electric held) — the canonical picture of the linear Stark 
ehect (LSE) splitting. The distance from the line center is expressed in a reduced scale (in 
the denominator there is the held strength in the corresponding ion-atom distance). The 
blue and red components are marked by proper colors, the central component - by black. 













Fig. 2.— Evolution of the intensity of selected Stark components of the line as a function 
of the reciprocal ion-atom distance in different approximations. For details — see text. The 
blue and red components are marked by proper colors. 
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Fig. 3.— Comparison of the Hq, line profiles of different references and meanings in the given 
physical conditions of plasma (rrip = 40 — the atomic mass of argon). For details, see text. 
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a;^(^ 

Fig. 4.— Comparison of the line profiles of different references and meanings in the given 
physical conditions of plasma. The symbols and the curves used have the same sense as in 
Fig. 131 The central dip displayed on the line-prohles represented by the dashed curves is too 
deep. This defect is the consequence of the negligence of the ion dynamics. The agreement 
of the measured and the calculated mFCSM data is excellent. 
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Fig. 5.— The pressure shift (PS) to the red of the peak of the Hq line as a function of 
the electron concentration in plasma. The dashed curve — the PS as in ApJ’87; the solid 
curve — PS of o ur present ca l culati ons using mFCSM; the points: the black triangles — 
measurements bv IWiese et al.l fjl972l) . the full circles — our measurements, both according 


to the paradigm of ETC. 
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Fig. 6.— The pressure shift (PS) to the red of the central dip of the line as a function of 
the electron concentration in plasma. The lines and the symbols have the following meanings: 
the black triangles and the dotted line — the original experi mental measurements of the ELC 
and the proper best £t, respectively, bv lWiese et al.l (119721) : the solid line and full circles — 
the ELC-type red shifts, calculated (via mFCSM) and measured, respectively, in the present 
paper; the short-long-dashed curve — the pressure red shift of the central dip according 
to the mFCSM calculations in the present paper; the short-dashed curve — the PS of the 
central dip as in ApJ’87. The points, the shapes of which are not commented on here, are 
taken from different bibliographical positions and, as a rule, belong to foreign authors. 


















Fig. 7.— The Stark shape (the dashed line) and the hnal effect of the Stark-Doppler 
convolution (the solid line) in one of the shallowest WD atmospheric layers. We see that 
the hnal shape, given by the convolution operation, differs completely from the initial Stark 
shape. This is of the Doppler-type one, without any central structure. 
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Fig. 8.— The Stark H/j shape (the solid line) and the hnal effect of the Stark-Doppler 
convolution (crosses) in one of the deepest WD atmospheric layers. We see that the shape, 
given by the convolution operation, is identical with the initial Stark shape. In deep, hot 
and dense atmospheric layers, where Stark broadening dominates totally, the contribution 
of the Doppler effect is completely negligible. 
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Fig. 9.— Superposition of the convolved Stark-Doppler profiles, corresponding to the suc¬ 
cessive WD atmospheric layers — from the outer most layers (the narrow, symmetrical, and 
unshifted shapes), down to the deepest ones (the extremely broad, asymmetrical, and red 
shifted shapes). The line-shape evolution is monstrous, indeed. 
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Fig. 10.— The same superposition as in Fig. [9l but, for clarity, containing only a few line 
profiles (H^ opacity shapes) from selected — every hfth — atmospheric layers. 
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Fig. 11.— The synthetic spectrum in the region of the line of a virtual WD of the given 
parameters, in the fluxes scale (left) and in the reduced (0-1) scale (right). 
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Fig. 12.— The Stark-induced shift (the pressure shift, PS) of the Balmer line, in ALC 
paradigm, of the WD atmosphere of the given parameters, as a function of the distance from 
the line center, in the fluxes scale. This hgure has been generated on the basis of the diagram 
of the type of Fig. ITTh . The solid line shows the present result, the thin dashed line — that 
in ApJ’87. For details, see the text. The dashed (horizontal) line represents approximately 
the leve l of the gravitational r ed shift in WDs of th e parameters, as in the legend above — 
cf., e.g. Provencal et ah (1998), Falcon et ah (2010). 
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Fig. 13.— The main results of the present paper — the Stark-induced red shifts, expressed in 
the Doppler velocity units, as functions of the distance from the line center, in the synthetic 
spectra of Hq, and Balmer lines of the WD of the given parameters. The diagram has 
been plotted on the basis of the synthetic Balmer line-prohles, expressed in the residual 
intensities, r^, in (0-1) scale of the type of Fig. ITTb . For detail, see the text. 






